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Abstract 

The effect of diffusional relaxation on the random sequential deposition pro- 
cess is studied in the limit of fast deposition. Expression for the coverage as a 
function of time are analytically derived for both the short-time and long-time 
regimes. These results are tested and compared with numerical simulations. 
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Diffusional relaxation in irreversible deposition processes of extended particles have been 
investigated recently In the deposition process of immobile particles, terms random 

sequential adsorbsion (RSA) the asymptotic behavior is dominated by the formation of 
'holes' too small for a new particle to fill in, resulting in the jamming of the available area. 
Diffusion allows non-effective depositions to be corrected in latter stages of the process. 
Thus, the coverage reaches its maximal value - the closest packing value - for long times. 
The asymptotic filling process is dominated by diffusive, power-law, approach to the steady 
state, as opposed to the exponential convergence to the jamming limit in the immobile case 
of lattice models. 

Mainly, two kinds of relaxation mechanisms were studied. Firstly, particle detachment 
has been studied, both experimentally and analytically. Exact results were obtained for 
special values of the parameters (detachment rate equal to deposition rate) ||. In many 
experimental situations, another relaxation process, diffusion of the deposited particles, is 
more significant. It has been shown numerically that the effect of this diffusional relaxation 
process in ID is a 1/y/i asymptotic approach to the closest packing value. This result 
was supported by analytical arguments as well [|J. A somewhat different model was also 
considered, in which the absorbed dimers are allowed to dissociate into two independent 
monomers. Each monomer can diffuse to its nearest neighbor sites. For this version of the 
model, for special values of the parameters (deposition rate twice the diffusion probability) 
exact solution is available |J. 

In this letter we study the combined effect of deposition and diffusion in ID, in the 
regime in which separation of time scales is possible - the deposition rate is high, and may 
be considered infinite with respect to the diffusion process. This regime is highly relevant 
for experimental interest. We find two series expansions, relevant for short and long times, 
which describe the whole dynamics of the filling process. 

Our model is defined as follows. We take an initially empty linear ID lattice containing 
TV sites, with periodic boundary conditions to minimize finite-size effects. First, particles 
are deposited randomly on the lattice up to the jamming limit. Each deposited particle fills 
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one lattice site, and excludes further deposition in his nearest neighbors sites (This is also 
equivalent to the deposition of dimers with no neighbor exclusion). Since the deposition is 
considered fast, this process takes zero time. Then, in each step one particle is selected, and 
with probability e/2 moves to the right or left, if possible. Each N steps are considered a 
time step, and thus the probability for a particle to move is e per time unit. Practically, 
when simulating the model numerically, we can select at each time step only eN particles 
and move them with probability 1/2 to the right or to the left After every movement, if 
a space for another particle was formed it is immediately deposited. We then look at the 
density as a function of time. 

In order to simplify the following analysis, we first set up our terminology. At any stage of 
the process, the lattice is filled by ordered regions, in which the particles are densely packed 
such that there is only one empty site between adjacent particles. In the border between 
these regions, there are two sucssesive empty sites separating the areas. Each region is 
termed "/c-mer" were k is the number of particles in the area. The initial concentration of 
the /c-mers can be easily calculated to be 



The only particles that move in any time are those which are at the edges of the /c-mers. 
As a result of these motions the /c-mers change their lengths. When a monomer moves, 
succesive three empty sites are obtained, a new particle is deposited, and the /c-mer and 
/c'-mer at two sides of the monomer become one big (k + k! + 2)-mer. 

Therefore, any change in the coverage results from monomer movement. Consequently, 
short times behavior is dominated firstly only by the monomer concentration. In latter 
times, monomers formed by the destruction of dimers contribute to the coverage as well, 
and thus the dimers, trimers, ... concentrations also play a role. In the short time regime, 
the dynamics of the /c-mers concentration c& is dominated by the transitions of a /c-mer to 
a (k ± l)-mer. The other process of a unification of a /c-mer, a monomer and a /c'-mer to a 
long (k + k' + 2)-mer can be neglected to first orders in time, since it generates only 4-mers 




(1) 
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or longer chains, and these do not contribute to the density up to fifth order. We thus have 
the following rate equations 

Ci = e(— 2ci + c 2 ) 

4 = e(c fc _i - 2c fc + Cfc+i), fc>l (2) 

with the initial conditions ([I]). Succesive approximations can be obtained by truncating the 
equation system after n equations, fixing c n+ i at its initial value. The result for c\ is then 
exact for n orders, and the coverage, which is given by 

Ap(t) = e f ci(t)dt (3) 

is exact for n + 1 orders. We thus get the first 4 approximations 

A«p(*) = cl{et) + (-2c? + c°)^ + (5c? - 4c° + c°)^£ + (-14c? + 14c° - 6c° + 4)^ + 0(t 5 ) 

(4) 

Figure 1 presents a comparison of the first 4 approximations with the results obtained from 
a numerical simulation using a. N = 256K lattice. Clearly, the fourth order expansion 
approximate the real curve up to et — 1. 

In order to study the long time behavior, we change our point of view. We term each two 
adjecent empty sites between successive /c-mers a "hole". Diffusion of particles at the edges 
of the monomers is equivalent to the diffusion of these holes ||. At the course of the hole 
diffusion, when two such holes are on adjacent sites, a space for a new particle is formed, 
and after its deposition the two holes annihilate. We thus see that our model is equivalent 
to the model of iV random walkers on a lattice which annihilate each other when joined. 
This model is well known and was used especially to describe the dynamics of chemical 
reactions of the type A + A — > inert [0-0 . In what follows we apply the analytical treatment 
developed for the reaction-diffusion problem to our model, and obtain an asymptotic series 
for the density. The two series, the asymptotic one and the previous short times expansion, 
describe the entire time regime very well, as can be shown by numerical results. 
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We wish to map our model to a standard model of annihilating random walks for which 
a rigorous result is known. We note that there are some differences between this standard 
model and the diffusing holes: (a) when a particle at the edge of a k— mer moves, the hole 
changes its position by two lattice sites. Thus, the diffusion constant is four times larger, 
(b) The distance between adjacent holes is always an odd number of sites. In particular, in 
the initial configuration the distances are odd. (c) The annihilation process occurs whenever 
the distance between the holes is one site and not when they are on the same site as in the 
standard model. Accordingly, in the initial state the minimal distance between successive 
k— mers is three sites. However, one expects these two differences to have no effect for long 
times for which the behavior is dominated by holes far from each other. This is confirmed 
by our numerical results. 

Thus we consider a model of random walks (RWs), originally distributed randomly on 
the lattice, with density 

p = p(0) = 2(l- Pr )=e- 2 , (5) 

where p r is the jamming limit density. The RWs moves with probability e two sites to right 
or left. When two such RWs join, they are annilihated. This model was solved exactly [|5]||, 
and we here follow the derivation given by Spouge [|J. Define fikin) to be the probability of 
the kth RW to be at a distance n from the origin 

A( n )=/(i-pr(^). (6) 

/3-(n) is then defined by 

oo 

P-(n) := 5 nfi + 2j2H k /3k(n) = 5 n>0 - 2p(l - 2p) n ~ l . (7) 
fc=i 

One also defines a(t; n) which is the probability that two RW whose original distance was n 
have meet until time t. In our model 

oo 

a(t;n) = I n (4Dt)+2 ]T h{Wt) (8) 

fc=n+l 
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where I n (x) := e~ x I n (x), I n (x) is the modified Bessel function of integer order, and D is the 
diffusion constant of the particles. 

Given the above definitions, Spouge's main result is 

pit) = pJ2a{t;n)f3^n (9) 

n 

Substituting the above expressions for our model we get 

p(t)/p=l-2pjt(l-2p) n [I n+1 (4Dt)+2 £ I k (ADt)} 

n=0 k=n+2 
oo 

= 1 - 2p J2 WDt) [q"' 1 + 2q k ' 2 + ... + 2g°] 

k=l 

oo 

= l-J2(2-q k - 1 -q k )I k (ADt) 

k=l 

oo 

= I Q (ADt) + (1 + q) £ q k -%{±Dt) (10) 

k=l 

where q = 1 — 2p. A similar expression was given by Balding et al 0. The particle density 
is given in terms of the hole density through the relation p = (1 — p)/2, and thus one obtains 
for the difference between the density and the maximal, closest packing, density 

Pc P - p{t) = p(t)/2 = — [h{Wt) + (1 + q) £ q'-'h^Dt)] (11) 

fc=i 

The diffusion constant is determined easily through the relation < r 2 >= 2Dt resulting 
in D = 2e. Now, the asymptotic behavior follows from the known asymptotics of the Bessel 
functions JL0|. 



i k {z) = - J=[l - ^ + (/i ~ .V * ~ 9) - ■■■], P = ^ (12) 



To the first order, all the Is are identical (k independent) and one has 

p M 1 + q, 1 0.0705.... , . 

In a similar way, next orders can be extracted. For example, the next correction is 

"-- p(t) = 8^3-(^ + ° ((rfr5/2) - (14) 
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where 



2 e 4 _ 4 e 2 + \ 

a 2 = — = 0.08886... (15) 



Figure 2 presents a comparison of the asymptotic leading order ( |T3| ) and the series ( [X0|) 
with numerical results obtained from a lattice of N = 1.2M sites. One sees that the whole 
series fits the results even for small values of et down to et = 0.1. 

In summary, it has been shown that the ID deposition-diffusion process leads to full 
coverage. The short time dynamics is determined by the temporal momnomer concentration. 
A fourth order expansion is given, valid up to et = 1. The long time kinetics is dominated 
by the attachment of two relatively long fc-mers which forms one long fc'-mer. This process 
is equivalent to the dynamics of a reaction-diffusion process, or to the probability of a RW 
to return to the origin. Thus, the asymptotic approach of the density to its saturated value 
in ID is 0(1/^/1). We derive an asymptotic series, based on this equivalency, valid for the 
intermediate and long time regimes (et > 0.1). 
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FIGURES 

FIG. 1. Numerical results for the coverage for short times, plotted vs. et, compared to the first 
four short-times approximations. 

FIG. 2. Numerical results for the coverage for long times, plotted vs. et, compared to the 
asymptotic leading order, and the Bessel functions' sum . 
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Density vs Time; Short time filling 
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Density vs. Time; Long time Filling 
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Eq- (13) 
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Eq. (11) 
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